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Abstract 

The origin of stabilization of complex microstructures along macrotwin 
boundaries in martensites is explained by comparing two models based on 
Ginzburg-Landau theory The first model incorporates a geometrically non- 
linear strain tensor to ensure that the Landau energy is invariant under rigid 
body rotations, while the second model uses a linearized strain tensor under 
the assumption that deformations and rotations are small. We show that the 
approximation in the second model does not always hold for martensites and 
that the experimental observations along macrotwin boundaries can only be 
reproduced by the geometrically nonlinear (exact) theory. 



1 Introduction 

Martensitic phase transformations are displacive first order transitions that involve 
a shear dominated change of shape in the underlying crystal lattice from a high 
symmetry phase, the austenite, to a degenerate low symmetry phase, the marten- 
site. Generally, the martensitic transition leads to very complex self-similar pat- 
terns that consist of twinned laminates on a wide range of length scales. It has 
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long been known that an important mechanism at the origin of these complex 
morphologies is the lattice mismatch between the different variants of the low- 
symmetry martensite and the austenite. Specifically, lattice continuity at the in- 
terfaces between two variants and between any variant and the austenite enforces 
long-range elastic interactions. The resulting strain energy can relax only through 
specific arrangements of twinned laminates that lead to a stress-free state charac- 
terized, at the scale of the laminates, by an invariant plane strain and, at a larger 
scale, by a very small average strain (see for example [1] for a general review). 

The dynamics of the transition is very often modeled through the use of a 
Time Dependant Ginzburg Landau theory which typically postulates that strain 
rates are proportional to a driving force derived from a thermodynamic poten- 
tial supplemented by compatibility-induced long-range elastic interactions. This 
overdamped description indeed generates complex microstructures, because it in- 
corporates a self-accomodation mechanism through highly anisotropic long-range 
elastic interactions. However, when the primary order parameter that controls the 
transition is an elastic long wavelength mode, as in proper martensites, this over- 
damped description is inadequate, as it does not incorporate the correct dynamics 
of the long- wave length sound waves. Indeed, using a simple linear analysis, it is 
easy to show that, because of the conservation law of density, the long-wavelength 
dynamics in the presence of damping leads to 

v 2 

where v s is the sound velocity, A some typical elastic constant and 7 a damping co- 
efficient. For sufficiently long wavelength and finite sound velocity, the period of 
oscillation l/v s q is always smaller than the lifetime \/v 2 s ^q 2 . The corresponding 
waves are underdamped and, therefore, propagate. In other words, we must con- 
sider that the sound velocity is finite and that strain-induced elastic interactions 
cannot reach their long-range character instantaneously. A correct description 
of the dynamics of proper martensites must therefore incorporate inertial effects, 
which requires the use of an inertial Lagrangian description. 

Another important feature of martensites is the presence of lattice rotations. 
More precisely, the transformation strain from austenite to martensite may be very 
large (up to the order of 10%) and thus lattice rotations may be important. This 
requires the use of a model that does not penalize rotations and hence is rotation- 
ally invariant. A model is considered to be invariant under rotations if its energy 
and dissipation potential are unchanged when its deformation field is rotated (see 
for example [2, 3] and the references therein). Most dynamic simulations of the 
martensitic transition have used strain-based models that are geometrically linear, 
i.e. that assume a linear relationship between the displacement and strain fields. 
This linearization simplifies greatly the formalism and allows the use of simple 



2 



and stable numerical schemes. However, linear geometry is obtained as an ap- 
proximation of the exact geometry under the assumption that displacements are 
infinitesimally small. As a consequence, the description is not invariant under 
rigid body rotations, in contrast to nonlinear geometry, which preserves the exact 
nonlinear relationship between the displacement and strain fields and is therefore 
rotationally invariant. 

Analytical studies of static energy minimizing geometrically nonlinear models 
and their geometrically linear counterparts indicate that there are significant dif- 
ferences between the two models and in idealized conditions the latter can even 
reproduce qualitatively incorrect results [1, 4, 5]. However, the effects of ge- 
ometrical nonlinearities on the dynamics in realistic conditions are still an open 
question. One may easily argue that neglecting these nonlinearities may have 
drastic consequences in an iterative process. 

In this paper, we present a numerical comparison of geometrically nonlinear 
and linear modeling for identical physical problems from thermodynamical, elas- 
tic and dissipative point of view and for the same physical parameters (mass den- 
sity, elastic constants, damping coefficient and interfacial energy). We show that 
dynamic models for pattern formation in martensites with geometrically nonlin- 
ear strains capture physically relevant features not captured in geometrically linear 
models. In particular, we show that the experimentally observed macrotwin inter- 
faces can only be reproduced within a nonlinear theory, as long-lived metastable 
states which cannot be reached in the geometrically linear model. 

2 Viscoelastic models for martensitic phase trans- 
formations 

A martensitic transformation is characterized by a displacement field u(x) = 
y(x) — x, where y{x) and x are the coordinates of a material point in the de- 
formed and undeformed states, respectively. Elastic energy F L depends only on 
the deformation gradient F = I + Vu and should be rotationally invariant. A sim- 
ple way to enforce this property is to assume that F L is only function of the La- 
grangian tensor e NL = (F T F — I) /2, which amounts to assuming, as usual, that 
the elastic energy depends only on the length changes between material points. 
This functional dependance guaranties that F L is invariant under rigid body ro- 
tations, simply because the strain tensor e NL itself is invariant if a rotation Q is 
applied to the deformed state: 

~e NL = \ ((QF) T (QF) - I) = \ (F T F - I) = e NL . (1) 



3 



We use the superscript NL to stress the fact that the Lagrangian strain is a nonlin- 
ear function of the displacement gradient : 

c nl _ 1 ( du i ! du j ! d^duk\ 
lj 2 \dxj dxi dxi dxj J ' 

The linear geometry approximation consists in neglecting the quadratic terms in 
the previous equation [1,6], leading to a linear strain tensor e L : 

i.e. e L = \ (F + F T ) — I. Obviously, e L is not rotationally invariant. 

In this paper, we consider the square to rectangle transformation as the sim- 
plest two dimensional model case for pattern formation in martensites (we will 
comment on that point below). The symmetry adapted linear combination of strain 
components are: 

ei e n + e 22 , e 2 e n - e 22 , e 3 e 12 e 2 \. (4) 

The quantity e 2 corresponds to the deviatoric strain and is the primary order pa- 
rameter, as it controls the square to rectangle transformation. The first and third 
terms e\ and e 3 correspond to the dilatational and shear strains, respectively, and 
play the role of secondary order parameters. To the lowest order, the correspond- 
ing strain energy density is given by the following Landau potential: 



F L = a 2 e\ + a±e\ + a§e\ + b x e\ + 6363. (5) 

The coefficients a 2 , a 4 , a 6 , 61, and 63 can be tuned to reproduce the elastic con- 
stants and the stress-free transformation strain of a specific alloy. We also include 
a gradient contribution, in the form of a Ginzburg term, in order to penalize in- 
terfaces between variants. This term of course breaks the scale invariance of the 
elastic medium and provides us with a length scale. For the sake of simplicity, we 
used here the following simple nonlocal density: 

F G = I (A«i) 2 (6) 

which fulfills the constraint of being rotationally invariant. Finally, we include a 
dissipation mechanism through a Rayleigh dissipation density given by: 

R = \^e- (7) 
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where, for simplicity, the viscosity coefficients {7^} are supposed to be indepen- 
dent of the amplitudes of the local order parameters e*. Taking into account a 
kinetic energy density T = \piiu we construct the Lagrangian density L = 
T — F L — F G . Finally, the dynamics of the system is given by the Langrange- 
Rayleigh equations: 

d dL dL OR 



dt diti dui diti 
These dynamical equations may be written in the general form: 



(8) 



„ _ 0F L _ d d dF G + 3 OR 

^ 1 dxj duij dxj dxk duijk 9xj duij 

where as usual u^j and Uij k stand for ^ and dx |^ , respectively. 

This formulation fits both for geometrically nonlinear and linear models. These 
two models differ only in the way the strain components enter, through the or- 
der parameters e*, into the definitions of the local strain energy F L and Rayleigh 
dissipation term R. They are linked to the displacement field u { by Eq. 2 for the 
nonlinear model and Eq. 3 for the linear one. 

We have chosen these equations as the simplest geometrically nonlinear mod- 
els, which can be compared to recent geometrically linear implementations of 
phase field (without inertia) [7, 8, 9] and Lagrangian dynamics [10, 11, 12] mod- 
els for martensitic phase transformations. In the following simulations, the dissi- 
pation term defined in Eq. 7 is treated linearly (i.e. with Eq. 3) for numerical sim- 
plicity. The effect of a nonlinear dissipation term is discussed elsewhere [13, 14]. 
However, we mention here that the models with linear and nonlinear dissipation 
terms lead to very similar final microstructures, simply because dissipation is van- 
ishingly small in the late stage of the dynamics. 



3 Experimental motivation 

We want here to investigate the dynamics and morphology of macrotwins, such 
as those observed in Ni 65 Al 3 5 [15, 16, 17], a shape-memory alloy which under- 
goes a cubic to tetragonal martensitic transition upon cooling (see Fig. 1). These 
configurations are very common in martensites, as they are one of the building 
blocks of the elastic energy minimizing microstructures generally observed in the 
late stage of the transformation. We will here specify the discussion to the cubic 
to tetragonal transformation, because this situation may be compared to the 2D 
square to rectangle transition modeled in this paper (see the discussion below). 

Generally speaking, when the martensitic transition proceeds, multiply twinned 
martensitic plates are formed to accomodate the shape change and minimize the 
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(a) 



(b) 



Figure 1 : (a) Typical TEM image of a macrotwin boundary in Ni 65 Al 35 revealing 
step type microstructures. (b) Typical HRTEM image of nanoscale martensite 
needles close to one of the steps of the macrotwin image shown in (a). Taken 
from [16]. 

elastic energy. Each plate consists of two of the three tetragonal Ll variants, 
separated by microtwin walls. The microtwin walls inside a plate originate from 
{110}-type planes of the B2 austenite. When the transformation proceeds, each 
plate continues to grow: a macrotwin plane forms whenever two such poly twinned 
plates come into contact and try to accommodate. We consider here the situation 
depicted schematically in Fig. 2, where two polytwinned plates are formed with 
alternating sequences of variants I and II, whose tetragonal axis are originally 
parallel to directions [100] and [010] of the austenite, respectively. A first plate 
is formed when variants I and II develop alternatively along microtwin walls that 
originate from former (110) planes of the austenite. As seen in Fig. 2(a), this ac- 
commodation mechanism requires rotations of variants I and II opposite in sign 1 
and which are not small if the shape change is large. More precisely, if the trans- 
formation strain U that transforms the austenite to variant I is written as \ 

E/ = Diag[)9,a,a] (10) 

the angle 9 is given by tan 9 = In the case of Ni 65 Al 3 5 [16, 17], the entries 

of the transformation strain have been estimated to a ~ 0.93 and /5 ~ 1.15, which 
leads to 9 ~ 6°. Another possibility for the same variants I and II is to accom- 
modate along twin walls originating from former (110) planes of the austenite, as 
in Fig. 2(b). This second configuration requires rotations opposite in sign to the 

during the growth process, compatibility of the plate with the austenite requires that there is 
a small rigid-body rotation of the entire plate that destroys this balance. 
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(a) 



(b) 



Figure 2: Twin walls involving variant I and II, whose tetragonal axis lie along 
directions [100] and [010], respectively, (a): the twin wall is parallel to a (110) 
plane of the austenite. (b): the twin wall is parallel to a (110) plane. In the case of 
Ni 65 Al 35 , we have 6 ~ 6°. 

ones involved in the first plate. 

Now, when these two microtwinned plates come into contact along a com- 
mon interface, this interface often reveals a zig-zag configuration consisting of 
successive well defined steps, as seen in Fig. 1(a). Along one of these steps, the 
microtwins of one plate end at a twin wall of the other plate (see Fig. 1(b)), this 
wall thus providing the macrotwin interface, which is locally a prior {110} plane 
of the austenite. At a larger scale, the steps alternate and we observe a zig-zag 
configuration of macrotwin interfaces parallel to prior (110) and (110) planes. An 
important feature of the macrotwin interface is the specific shape adopted by the 
twins that are perpendicular to the boundary: as seen in Fig. 1(b), alternate twins 
narrow and bend when approaching the macrotwin interface. 



Our aim is to analyze macrotwins associated to the cubic-to-tetragonal transfor- 
mation, but here we consider the simpler 2D square to rectangle transition. The 
latter may be interpreted as a 3D cubic to rectangle transition with the follow- 
ing restrictions. Consider a 3D crystal where the displacement field fulfills the 
following constraints: components u\ and u 2 are translationally invariant along 
axis x 3 and the third component u 3 varies linearly with x 3 , i.e. ui(xi, x 2 , x 3 ) = 



4 Numerical Results 
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u 1 (x 1 ,x 2 ), u 2 {xi,x 2 ,x 3 ) = u 2 {x 1 ,x 2 ) andu 3 (xi,x 2 ,X3) = \x 3 . The strain ten- 
sor 6ij then adopt the following restricted form: 

(en(x,y) ^ 12 (x,y) \ 
e 2 i(x,y) e 22 (x,y) (11) 
\ + £J 

where we used the nonlinear formulation. Obviously, if A = a — 1, where a 
is defined in Eq. 10, any laminate composed of tetragonal variants I and II, and 
hence with {110}-type twin walls, will be stress free. Therefore, from the mi- 
crostructural and energetic point of views, this laminate is exactly equivalent to a 
2D laminate consisting of the two rectangular variants of the square-to-rectangle 
transition with {ll}-type twin walls. In other words, we conclude that our 2D 
model should reproduce correctly the morphology of 3D multiple laminate struc- 
tures formed by tetragonal variants I and II and the corresponding macrotwins, 
such as those observed in the late stage of the martensitic transition. Of course, 
we cannot expect that a 2D model will capture correctly the dynamics itself, in 
particular when austenite and martensite coexist. However, the results reported 
below concern specifically the stability of macrotwins in the late stage of the tran- 
sition. 

The evolution equations were solved using a Fast Fourier Transform method 
with implicit-explicit finite difference timestepping schemes [13, 18, 19]. Equa- 
tions are first written in a dimensionless form using adapted units for mass density 
(p = 6.657 g/cm 3 ), time (£ = 10" 13 s), spatial coordinate (d = 0.1 nm) and en- 
ergy density (f = 29.8 GPa). The dimensionless parameters used here are the fol- 
lowing: p=l 9 a 2 = -0.40, a 4 = -2.18, a 6 = 119, b x = 2, b s = 2, 7 = 0.1, j3 = 
0.1, where a refers to the dimensionless counterpart of the physical quantity a. 
These parameters have been chosen to make our 2D system representative of the 
situation in Ni 65 Al 3 5. This requires the knowledge of experimental quantities, 
such as the twin interfacial energy and, as we are concerned with microstructures 
in the martensitic state, elastic constants of the martensite. Unfortunately, as far as 
we know, these constants are not known. However, experimental measurements of 
elastic constants in the austenite are available. In particular, using ultrasonic tech- 
nics [20], it has been observed that the elastic constants C L = (C11+C12+2C44) /2 
and C44 varies only slightly when approaching the martensitic transition, contrary 
to the shear elastic constant C = (Cn — Ci 2 )/2. We thus assumed that the mea- 
surements for C L and C 44 , with C L ~ 290 GPa and C 44 = 120 GPa, can be used 
to estimate the corresponding elastic constants in the martensitic phase 2 . In addi- 
tion, as we consider here a situation which is far below the transition temperature, 

2 It is understood here that the elastic constants of the martensite are expressed in the unde- 
formed reference state, i.e. in the coordinate axis of the austenite. With the simple strain energy 
used here (Eq. 5), it is easy to show that, within this coordinate system, the elastic constants of 
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the shear elastic constant C should not be particularly small. We assume that 
it is of the same order as C 44 . Specifically, we used C ~ C44 (similar results 
have been obtained with C ~ 2C44). An estimation of the twin interfacial energy 
is also needed, in order to fix the length scale of the simulation. We chose here 
a ~ 75 mJm -2 . We mention that a specific estimate of a is needed only to fix the 
length scale of the system. Any choice leads to the same dimensionless equations 
and, therefore, to the same microstructural evolution. Finally, we mention that the 
choice presented above leads to an equilibrium stress free deviatoric strain close 
to the experimental one, namely e^ q = 0.2. 

The simulations are carried out on a grid of size 1024 x 1024 that corresponds 
to a physical size of 0.1/xmx0.1/im. The interfaces within the simulation spread 
typically over three grid points. This corresponds to an interface width of the order 
of 0.3 nm, which is close to typical interface widths observed in the experiments. 

Motivated by the experimental situation discussed above, which reveals struc- 
tural features corresponding to large crystal lattice rotations along macrotwin 
boundaries [16], the simulations start with the initial configuration shown in Fig. 
3. It consists of four macrotwins between laminates formed by alternating se- 
quences of variant I and II. The deviatoric order parameter field e 2 is initialized to 
small values, proportional to the equilibrium value of the corresponding variant, 
and e\ and e 3 are set identically to zero. The initial velocity of the displacement 
fields Ui is also set to zero. 




Figure 3: Initial condition for the deviatoric strain e 2 . The order parameter e 2 
is initialized in the form of sinusoidal waves of amplitude 0.01 (the equilibrium 
deviatoric strain is e 2 = ±0.2) and e\ and e 3 are set identically to zero. Red and 
blue domains correspond to variant I and II, respectively. 

The dynamics is first run within the exact nonlinear geometry scheme. The 

the martensite verify C\\ = C 22 . An inequality would be obtained if a coupling between the 
dilatational e\ and deviatoric e 2 strains was included in the strain energy density. 
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(a) (b) (c) 



Figure 4: Microstructure evolution in the geometrically nonlinear model. The 
reduced times from (a)-(c) are t = 40, 100, 400 x 10 2 . 

corresponding microstructural evolution is shown Fig.4. We observe that the mi- 
crotwins in the neighborhood of macrotwin boundaries display new features at 
early stages. In particular, a variant that meets a perpendicular lamella made of 
the same orientational variant is seen to bend and to narrow when approaching the 
macrotwin interface (see Fig. 4(a)). Finally, as seen in Fig. 4(c), the final state 
obtained using the nonlinear modeling still displays the initial macrotwin bound- 
aries, but characterized now, on one side, by alternating sequences of variants that 
narrow and bend when approaching the macrotwin boundary. We mention that 
this final state, although metastable, does not evolve anymore. 

For comparison, we present in Fig. 5 the microstructural evolution within 
the linear geometry. We used the same material parameters as for the nonlinear 
modeling. The only difference is that we used here Eq. 3 instead of Eq. 2. We 
note that the final microstructure is completely different from the one obtained 
with the linear model. We will comment further on this point below. 

In Fig. 6, enlargement of the circled region on Fig. 4(c) is shown together 
with a HRTEM image of a macrotwin boundary in NiAl [17] for a convincing 
comparison between the nonlinear modeling and the experimental observation. 
We indeed observe that the tapering and bending of alternating variants is well 
reproduced by the nonlinear model. We note also that needle splitting and hook- 
type forms are also reproduced. 

5 Discussion 

The driving force for the development of these microstructures is of course the 
reduction of the strain energy. The underlying mechanism is the following. Con- 
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Figure 6: Comparison between simulated microstructure and experiments, (a) 
Enlargement of circled region in Fig. 4(c), rotated by 45°. (b) HRTEM image of 
a macrotwin boundary [17]. 
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Figure 7: Schematic representation of a macrotwin boundary between two differ- 
ently oriented laminates. 

sider a macrotwin interface parallel to a (110) plane, as the one sketched in Fig. 
7. This situation corresponds to the macrotwin boundary in the upper right corner 
of the final configuration displayed in Fig. 4(c) and also in Fig. 6(a). Accom- 
modation, along the (110) twin plane, of variant I (on the right side of the plane) 
and of variant II (on the left side) require rotations equal to —6 and +0, respec- 
tively. Now, as variant II has already been rotated by +0, accommodation along 
the former (110) twin plane that collides with the macrotwin boundary requires 
that variant I rotates by 30 in order to fit coherently with variant II. As a result, 
the former (110) twin wall must rotate by 20. This is at the origin of the bending 
of the twin walls that we observe in the simulation as well as in the experiment. 
We note also that the different rotations displayed by the variants as well as the 
amplitude of the bending observed in the simulation are numerically close to the 
expected values (see Figs. 6 and 8(b)). Another important feature of the previ- 
ous accommodation mechanism is that it leads, along the macrotwin boundary, to 
interfaces between domains made of the same orientational variant (i.e. domains 
with identical deviatoric strain e 2 ) but with very different rotations. Specifically, 
in the situation depicted in Fig. 7, there is an interface between two domains of 
variant I that have been rotated by —6 and 30, respectively, leading to a high strain 
energy density along this interface. The system will try to shorten this high energy 
interface and, consequently, to lengthen the neighboring low energy twin walls. 
This mechanism is at the origin of tapering of alternate variants observed in the 
simulation and in the experiment shown in Fig. 6. 
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In brief, the bending and narrowing of alternate variants along a macrotwin 
boundary requires large local rigid body rotations (as high as 30, with 6 ~ 5°). 
This illustrates the importance of using a rotationally invariant modeling, i.e. a 
geometrically nonlinear scheme. 

Indeed, as shown in Fig. 5, where we report the microstructural evolution 
within the linear geometry, the time evolution with linear model is completely 
different. As in the nonlinear modeling, stress accommodation operates first at 
short length scales, as can be seen in Fig. 5(a) where we observe the formation 
of needles in the neighborhood of the macrotwin boundaries. However, in further 
time steps, the needles retract and martensite laminates start to disappear in favor 
of large martensitic domains (Fig. 5(b)). Finally, the geometrically linear model 
produces a final state that extends across the whole domain, i.e., a simple laminate, 
shown in Fig. 5(c). 

The differences between the microstructures obtained within the two models 
can be directly linked to the role of rotations involved by the accommodation 
mechanism. For a quantitative understanding of this point, we present in Fig. 8 
the variation of the strain energy density of an homogeneous variant as function 
of the angle (f) of an applied rigid body rotation, according to the linear model, i.e. 
using Eqs. 4 and 5 with the linear relation of Eq. 3. We remind that this strain 
energy density does not depend on (j) if a nonlinear scheme is used. The situation 
is of course very different within the linear model. The strain components e^, 
and therefore the order parameters e*, are no longer rotationally invariant. In par- 
ticular, a rotation of angle (f) will generate non-zero dilatational and shear strains 
ei and e 3 . As seen in Fig. 8, the strain energy density difference between the 
martensite and the austenite, changes very slowly with (p for small angles, say 
up to approximately 10°. However, most of the stabilization of the martensite is 
lost for higher rotations. Indeed, rotations of the order of 30, as those involved in 
the stabilization of the needles within the nonlinear model, almost exhaust the en- 
ergy difference between martensite and austenite. This unphysical behavior is the 
reason why the dynamics within the linear model cannot sustain the bending and 
tapering of the microtwins and, therefore, cannot stabilize the macrotwin bound- 
aries observed in the nonlinear model and in the experiment discussed above. 

Another convincing illustration of this point is presented in Fig. 9, where we 
show the strain energy density maps of the final state obtained within the exact 
(nonlinear) model. Two maps are presented: (a) the strain energy density for the 
final (metastable) configuration obtained with the nonlinear model and, (b), the 
strain energy density for the same configuration, but calculated using the linear 
model, i.e. by replacing Eq. 2 by Eq. 3. The local strain energy concentration 
in the exact model is locally concentrated on variant tips in very small regions 
(Fig. 9(a)), precisely where the needles taper against the macrotwins, leaving short 
but highly energetic interfaces between identical variants. However, this high 
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Figure 8: (a) Strain energy of an homogeneous variant as a function of the angle 
of a rigid body rotation in the linear model, (b) Simulated microstructure at grid 
scale in the nonlinear model. This snapshot is the magnification of the macrotwin 
boundary shown in Fig. 6(a). The angle (j) measures the rotation of variants with 
respect to the reference state. 



energy concentration is not strong enough to create a driving force to move the 
domain boundaries and consequently this metastable state can persist in the non- 
linear model. In contrast, energy concentration is high on very large regions in 
the linear model (Fig. 9(b)), precisely where finite rotations are needed to accom- 
modate variants along twin walls. As a result, the excess strain energy in the 
intermediate layer, where alternate variants bend, are too large for bending and 
splitting type microstructures to be long lived metastable states and therefore they 
do not persist in the linear model. 

There are other computational works in the literature that claimed to observe 
similar microstuctures to those obtained here using the nonlinear model but none 
of them compared geometrically nonlinear and linear models [7, 21, 22] . For 
example, in [21], a geometrically nonlinear Ginzburg-Landau approach was used 
to analyze twin morphologies. Using a simple conjugate-gradient minimization 
algorithm, twin narrowing was observed at the junction of two mutually perpen- 
dicular monodomains made of the same orientational variant. However, large 
macrotwin boundaries were not analyzed and comparison with a geometrically 
linear model was not undertaken. Needle shape morphologies have also been 
observed in [7], using a simple geometrically linear dissipative Time Dependant 
Ginzburg-Landau approach. Under applied stress, nucleation of twins between 
previously existing variants was observed. This new twin generation exhibits nee- 
dle shapes, but before they collide with the pre-existing variants. Moreover, these 
microstructures were only short-lived and stable macrotwin boundaries between 
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Figure 9: (a) Strain energy map of the final microstructure obtained in the non- 
linear modeling. Black corresponds to the higher energy concentration, (b) strain 
energy map of the final microstructure obtained in the nonlinear modeling, but 
computed after linearization of the relation between displacement and strain. 



perpendicular laminates was not observed. [22] observed complicated long lived 
microstructures in a geometrically linear model, but no well-defined macrotwin 
interfaces as those analyzed here were found, and stability of multiple laminates 
was observed only if a specific defect field was introduced. Concerning that point, 
we note that, according to the experimental observations, macrotwin interfaces do 
not seem to be stabilized by defects [15]. 

6 Conclusion 

In summary, we have presented here a numerical comparison between two dynam- 
ical models for understanding the formation of complex microstructures observed 
in materials undergoing a martensitic phase transformation. The first model in- 
corporates a geometrically nonlinear strain tensor, to insure that the strain energy 
is rotationally invariant, whereas the second one is based on its geometrically 
linear approximation. The linear model is commonly used as an approximate 
approach for studying the dynamics of martensitic phase transformations. Our 
findings show that this approximation cannot capture and reproduce the physical 
mechanisms at the root of the stabilization of the late stage microstructures ob- 
served in martensites. In particular, we show that the model with geometrically 
nonlinear strains produces final metastable states with multiple laminates, sepa- 
rated by stable macrotwin interfaces along which bending, tapering and splitting 
of microtwin needles is observed, in agreement with experiment. On the other 
hand, the linear model produces final states that are always simple laminates. We 
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argue that the absence of stable macrotwin boundaries in the linear model is due 
to the penalty resulting from the energetically costly dilatational and shear strains 
generated, within a linear geometry scheme, by the lattice rotations required for 
the coherent accommodation of the deviatoric strain along macrotwins. 
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